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Directionality in the WMAP Polarization Data 
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ABSTRACT 

Polarization is the next frontier of CMB analysis, but its signal is dominated over 
much of the sky by foregrounds which must be carefully removed. To determine the 
efficacy of this cleaning it is necessary to have sensitive tests for residual foreground 
contamination in polarization sky maps. The dominant Galactic foregrounds introduce 
a large-scale anisotropy on to the sky, so it makes sense to use a statistic sensitive to 
overall directionality for this purpose. Here we adapt the rapidly computable T> statis- 
tic of Bunn and Scott to polarization data, and demonstrate its utility as a foreground 
monitor by applying it to the low resolution WMAP 3-yr sky maps. With a thorough 
simulation of the maps' noise properties, we find no evidence for contamination in the 
foreground cleaned sky maps. 
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servations - cosmology: theory - large-scale structure 



1 INTRODUCTION 

Our progress in understanding the large scale make-up of 
the Universe has been propelled rapidly in the last 15 
years by extremely detailed measurements of the Cosmic 
Microwave Background (CMB). Current measurements of 
the CMB temperature anisotropies have allowed us to place 
stringent constraints on some of our Universe's key pa- 
rameters and are proving to be a driving for ce in fun- 
damental physics (see e.g. IScott fc"s"m oot 2006). Cosmic- 
variance places limits on the amount of information which 
can be extracted from the temperature anisotropies alone, 
however. The Planck satellite, for example, is expected 
to reach this limit over the entire primary tempe rature 
anisotropy power spectrum iPlanck Bluebookl [2006). For- 
tunately, the polarization of the CMB anisotropies pro- 
vides us with a wealth of additional cosmological infor- 
mation which we have only just begun to probe. The 
benefits of CMB polarization data are numerous. Perhaps 
most excitingly, polarization data provide a chance to de- 
tect primordial gravitational w aves, through the measure- 
ment of B-mode polarization ( Sgljak & Zaldarriaga 19971: 
Kamionkowski, Kosowskv fc Stebbins I Il997l ; Hu fc White I 



19971 : iKamionkowski fc Kosowskv! 1 199S ) . The rapid im 



provements in CM B polarization data provide d by experi- 
ment s like PIQUE (iHedman et al.l200ll). D ASI i|Kovac et all 
120021) CAPMAP llBarkats et all |2005|), B OOMERanG 
Masi et alj|2006l ), CBI (Rcadhcad et al. 2004), MAXIPOL 



Wu et al.ll2006l ). WMAP |Page et al. 2006), and upcoming 
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experiments like BICEP, CLOVER (|Tavlorll2006l ), QUAD, 
and the Planck satellite, ensure that polarization will be- 
come one of the most important tools we have for furthering 
our understanding of the Universe in the years to come. 

Unfortunately, the polarization signal is much more dif- 
ficult to measure than the temperature anisotropies. The 
magnitude of the CMB polarization anisotropies is at most 
10 per cent that of the temperature anisotropies, making 
their detection alone a challenging task. This difficulty is 
being gradually overcome with improvements in detector 
technology. Another problem with the CMB polarization, 
however, is that its signal is dominated over most of the 
sky by foregrounds. In particular, over the wavelength range 
where the CMB is typically measured, synchrotron and non- 
spherical dust emission from within our Galaxy contaminate 
uncleaned polarization maps. It is not yet clear how well the 
problem of foreground contamination can be solved. Its reso- 
lution will require both a detailed knowledge of the structure 
and magnitude of the foregrounds, and an understanding of 
how to remove them from the polarization data without in- 
troducing unwanted systematics. 

Understandably, the issue of polarization foregrounds 
has received much attention in the literature. Several pa- 
pers have been written on the construction of templates 
for the emission ide Oliveira-Costa et al] [20031 : iBaccigalupil 
l2003l : |Page et al.ll2006l : lHansen et alj|200q ). and a number of 
methods have been proposed for performing the foreground 
removal dBouchet. Prunet fc Sethil 19991 : lHansen et al ] |2006l : 
ISlosar. Seliak fc Makarovll2004l ). Testing the efficacy of these 
methods is a crucial step in the analysis pipeline, to compare 
the abilities of the various removal techniques and ultimately 
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to provide confidence in the accuracy of CMB polarization 
data. Tests for the existence of residual polarized foreground 
contamination have received less attention in the literature, 
however. The test used in the WMAP team's analysis, where 
foregrounds are cleaned using a template meth od in pixel 
space, is a simple \ 2 statistic l|Page et al.ll2006l ). Unfortu- 
nately, this statistic cannot be used to make comparisons 
between different cleaning methods. The only currently pro- 
posed 'general' statistic which is capable of this application 
is based on the B ipolar Power Spectrum (BiPS) and has 
been developed bv lBasak. Haiian fc Souradeep I (|2006l ). 

In this paper, we introduce and analyse another general 
method of testing for residual foreground co ntamination in 
CMB polarization maps, the T> statistic of iBunn fc Scottl 
( 2000) . It possesses the virtues of having a simple interpreta- 
tion as a measure of directionality in a CMB map and of be- 
ing extremely rapid to compute. We argue that this statistic 
is well suited to the detection of foreground contamination, 
as both the synchrotron and dust foreground components 
have strong overall directionality on large scales. We then 
apply the T> statistic to the WMAP 3-yr polarization maps 
to demonstrate its usefulness for real data. 



2 THE V STATISTIC 

The T> statistic was first presented bv lBunn fc Scotd l|2000h : 
the V value for a sky map pixelized into N pixels is defined 

as 

p= max ft /(n) 
mnift / (n) 

where n runs over the unit sphere, and with the function 
/(n) given by 
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/(«) = ^2 w p(n-g p ) 2 . 



(2) 



Here w p are weights chosen to compensate for noise and the 
effects of a cut sky, and g p is a local directionality vector 
calculated for the pixel p. The f(n) value is a measure of 
how much (or how little) the g p of a map tend to point 
in a given direction. The T> statistic tells us how extreme 
these values are. Once T> has been calculated for real sky 
data, we can compare its value to that found for simulated 
realizations of the CMB. Calculating T> for a large number 
of these simulations gives us a histogram of expected values, 
and excess directionality in a CMB data-set appears as a 
value of T> which is an outlier of this distribution. 

For sky maps of temperature, a reasonable quantity to 
use for g p is the local temperature gradient VT P . For polar- 
ization maps, on the other hand, the polarization psuedo- 
vectors themselves are obvious candidates. We define g p to 
be a vector whose magnitude and direction are those of the 
polarization at point p. To express g p explicitly in terms of 
the Stokes parameters Q and U, recall that the polarization 
magnitude is 



VQ 2 + u 2 , 



(3) 



while the polarization direction lies in the tangent plane 
to the celestial sphere at p and makes an angle with the 
meridian of 



Measurin g the angle from the meridian is a WMAP con- 
vention (jPage et al.ll2006l ) which we adopt here, although it 
is not universally followed. We define the vector g p to have 
magnitude P and direction given by 6. 

Because the polarization is a spin-2 quantity, it is repre- 
sented by headless 'pseudo- vectors' rather than vectors, so 9 
can always be rotated by 180° without altering the polariza- 
tion. However, due to the quadratic definition of /(n), this 
ambiguity does not affect T>. Moreover, because the func- 
tion / smoothly averages together polarization information 
over the whole sky, T> is only sensitive to large scale (> 20°) 
directionality. 

The weights w p must be chosen to satisfy the null hy- 
pothesis that for isotropically distributed g p no preferred 
direction will be found on average, even when parts of the 
sky have been removed. We can see how these weights may 
be obtained by casting f(n) in matrix form as 

f(n)=n T An, (5) 

where A is a 3 x 3 matrix defined by 



Ai 



N 

— S w p g p ig p j. 
P =i 



(6) 



In order for no preferred direction to be found, we must 
have f(n) independent of n (on average). From equation[5]it 
can be seen that this constraint may be realized by choosing 
suitably normalized w p such that (A) is equal to the iden- 
tity matrix for isotropically distributed g p . This criterion 
ensures that /(n) will be approximately constant over the 
sky when the underlying data lack a preferred direction, al- 
though it does not guarantee that fluctuations of / about its 
ensemble average will be isotropic. As a result, D may still 
pick out preferred directions arising from the noise structure 
or sky coverage, even when the underlying signal is isotropic. 
We find that such departures from anisotropy are weak in 
practice; this point is discussed further in Section [3] 

Because A is symmetric, requiring (^4) = I gives only 
6 independent equations on the weights. To completely de- 
termine them for a large number of pixels, the additional 
constraint is imposed that Va.r(w p P p ) be minimized, where 
P p is the expected value of \g p \ 2 given by pixel noise and 
cosmological signal. For a large signal-to-noise ratio on an 
isotropic sky this is equivalent to minimizing the variance 
of the weights. When noise is not negligible, minimizing 
Var(w p P p ) rather than Var(ui p ) has the favourable effect 
of down- weighting noisy pixels. More explicit details on th e 
calculation of the weights are given in lBunn fc Scottl (|2000h . 

Once the weights w p have been calculated, the simple 
quadratic definition of f(n) makes it possible to calculate 
T> for a given data-set in O(N) operations, where N is the 
number of pixels. This can be done by noting that the values 
of n which extremize /(n) are given by the eigenvectors of 
A. T> follows immediately by taking the ratio of the largest 
and smallest eigenvalues. 

D is one of the simplest statistics which can be pro- 
posed for the identification of statistical anisotropy in 

a CMB map. To compare it to the BiPS method of 

Basa kT"rIauan fc Souradeep I (|2006l ). we note that it is ul- 
timately a less powerful statistic, as BiPS can be used to 
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probe both large and small scales, whereas T> is only sensi- 
tive to large scale anisotropies. The advantage of T>, however, 
is that it can be calculated in O(N) operations, an unbeat- 
able scaling. This speed makes T> an excellent statistic for 
'first look' determination of anisotropy in CMB maps, and 
the concomitant possible foreground contamination. 



3 DIRECTIONALITY OF POLARIZATION 
FOREGROUNDS 

The two main sources of foreground contamination in CMB 
polarization maps at the ~ 100 GHz frequencies where the 
CMB is usually probed are Galactic synchrotron and ther- 
mal dust emission. Polarized synchrotron emission results 
from the acceleration of cosmic-ray electrons as they orbit in 
the Galactic magnetic field. Polarized dust emission, on the 
other hand, is the result of non-spherical dust grains which 
tend to align their long axes perpendicular to the field; these 
dust grains preferenti ally emit thermal radiatio n polarized 
along their long axis (|Davis fc Greensteinlll95lT ). 

The polarization of both the synchrotron and thermal 
dust emission has its origins in the anisotropy introduced by 
the Galactic magnetic field, and at CMB wavelengths (where 
Galactic Faraday rotation is negligible) these two compo- 
nents should both be polarized preferentially in the same 
direction. Thus, we expect that they have a common and 
distinctive directionality on the sky, dictated by the large 
scale structure of the Galactic magnetic field. To illustrate 
the signature of this directionality under application of the 
T> statistic, we analyse the low resolution K-band polariza- 
tion map at 23 GHz and the dust polarization template of 
the WMAP team, masking both with the standard polar- 
ization mask (P06) to mimic the situation when analysing 
maps for foreground contamination. For the K-band map 
we find a preferred direction at a Galactic latitude 86.4°, 
within 4° of the poles. For the dust emission template we 
find that the preferred direction lies at 79.1°, 11° shy of the 
poles. Thus, the significant presence of synchrotron and dust 
foregrounds in a polarization map is expected to result in an 
outlying value of T>, with a preferred direction in the vicinity 
of the Galactic poles. 



4 ANALYSIS OF WMAP 3-YR POLARIZATION 
DATA 

The release of the WMAP 3-year data has provided us with 
a glimpse at the structure of CMB polarization. Here we 
analyse the Q and U polarization maps for foreground con- 
tamination using the T> statistic. 

We have chosen to examine the WMAP low resolution 
3rd year (v2) polarization mapi|3 at HEALPl>jf] resolution 
4, corresponding to 3072 pixels on the entire sphere. These 
are the maps used for low i analysis of the WMAP data, cor- 
responding to the angular scales which T> effectively probes. 
At this resolution it is also feasible to accurately simulate 
realizations of the WMAP noise structure using the full 



1 Available at http://lambda.gsfc.nasa.gov/ 
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1.1 1.2 1.3 

Value of D 

Figure 1. Directionality histograms for maps simulated with 
WMAP low resolution V-band noise levels, with (solid) and with- 
out (dot-dashed) cosmological signal added. The two are essen- 
tially indistinguishable. The dashed line shown at T> = 1.37 is the 
value of T> calculated for the V-band data themselves. 



noise covariance matrix, which becomes necessary for large 
angle analyses. 

In all of our simulations, we generate cosmological 
signals given by th e ' WMAP +all ACDM' fit parameters 
|Spergel et al.| [2006). The details of this choice have little 
effect on the outcome of our simulations. All of the polar- 
ization maps are strongly noise-dominated, and we find that 
histograms of T> for simulated data with and without cos- 
mological signal are indistinguishable for all of the WMAP 
bands; this is illustrated in Fig. [T] for simulated V-band 
maps. Thus, we believe that, with an accurate treatment of 
the maps' noise properties, any excess directionality which 
we detect in the WMAP data can be attributed to residual 
foregrounds, and not a cosmological signal. 

For all of our simulations, we add correlated Gaussian 
noise to each unmasked map pixel based on the full QU co- 
variance matrices for the corresponding WMAP frequency 
band, as provided by the WMAP team. We generate this 
noise using the Cholesky decomposition of the covariance 
matrix. This method accounts for QQ, QU and UU corre- 
lations both within each pixel and with other pixels across 
the sky. 

For uncleaned maps, the T> statistic produces the ex- 
pected results. For each band, we generate and analyze 
~ 2000 simulated skies. An example of this is shown in Fig.rU 
We can use the number of standard deviations (at which the 
real data lie from the histogram mean) as an approximate 
measure of the significance with which we detect foreground 
contamination. This is illustrated in Fig. [2] which shows the 
expected dependence on radiometer band as the synchrotron 
emission decreases with frequency and the dust contribution 
increases, leading to an apparent minimum of foreground 
contamination near the V-band. We can see that use of the 
standard P06 polar mask results in a greater level of con- 
tamination than when we use the more stringent P02 mask. 
For the P06 mask, all of the preferred directions which are 
found for the real data lie within 10° of the Galactic poles, 
and for the P02 mask the preferred directions all lie within 
25° of the poles, further demonstrating the foreground origin 
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Figure 2. Significances of foreground contamination for the 
WMAP low resolution maps using noise based on the full co- 
variance matrix. The (□) symbols are for the standard P06 mask 
(27 per cent sky cut), while the (A) symbols are for the more 
stringent P02 masking (50 per cent sky cut for our downgrad- 
ing scheme). The curve roughly follows that expected from the 
frequency dependence of contamination by synchrotron and dust 
foreground components. 
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Figure 3. Directionality histograms for simulated realizations of 
the WMAP V-band (solid) and Q-band (dashed) 'low resolution 
foreground reduced' maps, using the full covariance matrix de- 
scription of the maps' noise properties. The vertical lines are the 
values of T> calculated for the Q- and V-band data themselves. 



of the directionality excess. With P06 masking, we find that 
98% of our V-band simulations have values of T> less than 
the value calculated for the WMAP data. With P02 mask- 
ing, this result drops to 90%. The corresponding significance 
figures of merit are 2.5<r and 1.3<r respectively. 

In order to assess the efficacy of the foreground removal 
performed by the WMAP team we also analyse the low reso- 



lution Q- and V-band foreground reduced maps. These have 
had dust and synchrotron templates projected out to remove 
foreground contamination. The results of these analyses are 
shown in Fig. [3] From this figure it can be seen that the 
Q and V maps have similar noise properties under analy- 
sis with T>, as the simulation histograms almost completely 
overlap. We find no evidence for residual foreground con- 
tamination in these maps, their T> statistics both lying well 
within the range of probable values, and with the maximal 
T> directions for the Q- and V-band data being 55° and 27° 
from the pole, respectively. Our significance figure of merit 
is 0.6cr for the Q-band and 0.16cr for the V-band maps. 31% 
of the simulated T> values lie below the map value for the Q- 
band, and 60% of the simulated data lie below the V-band 
value. Thus, we find no evidence for foreground contamina- 
tion in these cleaned maps. 



It is interesting to note that if we had generated noise 
which only accounted for the intra-pixel or diagonal compo- 
nents of the noise covariance matrix then we would have 
found strong evidence for residual contamination. In the 
case of the V-band foreground reduced data, for example, 
we would have found some evidence for foreground contam- 
ination, with over 98% of simulated T> values falling below 
that calculated for real data in the cleaned V-band maps 
(corresponding to a value of 2.5<r for our significance fig- 
ure). Accounting for inter- pixel noise correlations broadens 
the distribution of simulated T> values. Indeed, the detailed 
noise structure of the WMAP data appears to be such that 
it introduces a degree of anisotropy into the maps. This can 
best be seen by looking at the histograms of preferred direc- 
tions found in a series of 30,000 V-band map simulations, 
given in Fig. 2] Differing methods of noise generation were 
utilized for three sets of simulations. 'Diagonal' noise was 
generated from only the QQ and UU elements of the covari- 
ance matrix, intra-pixel noise included QU terms, and the 
full covariance matrix treatment simulated inter-pixel corre- 
lations as well. Unlike the isotropic distribution of preferred 
directions which is found when accounting only for diagonal 
noise correlations, there is a tendency towards anisotropy 
for the intra-pixel and full covariance simulations, with a 
definite slight preference for directions in the vicinity of 
(l,b) — (150°, 40°) for the full covariance simulations. The 
statistics of these two dimensional histograms support this 
visual finding, with the diagonal, intrapixel, and full co- 
variance simulation methods having variance/mean values 
of 1.41, 1.66 and 3.80 respectively. For an isotropic sky, we 
expect the histogram values to follow a multinomial distribu- 
tion, and to have a variance/mean of 1±.08. We believe that 
the excess variance which is seen in the diagonal and intra- 
pixel simulation histograms can be attributed to the slight 
inherent bias possessed by T>, as discussed in Section [2] The 
much higher excess variance seen in the inter-pixel simula- 
tions can be attributed to the sum of this bias and another, 
stronger bias given by the noise covariance structure of the 
map. This bias is independent of the specifics of the T> statis- 
tic, and most likely can be detected in the future by other 
polarization statistics sensitive to anisotropy. The depen- 
dence of our results on the precise treatment of noise un- 
derscores the general need for care when dealing with large 
angle polarization data. 
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Figure 4. Histograms of preferred direction found for maps simu- 
lated with diagonal, intra-pixel, and full covariance matrix based 
noise (Mollweide projection, Galactic coordinates). The colour 
units are in number of simulations. The preferred directions cal- 
culated by T> are headless, and so we have taken the convention 
of choosing them to be in the Northern Hemisphere. 



5 CONCLUSION 



We have investigated the V statistic of lBunn fc Scottl (|2000l ) 
in the context of CMB polarization data, finding that it is an 
excellent general tool for quickly assessing the magnitude of 
foreground contamination in polarization maps. Under anal- 
ysis with T>, imperfectly removed foregrounds leave a tell- 
tale signature of excess directionality toward the Galactic 
poles. An analysis of the WMAP 3-year polarization data 
supports these statements. We find the expected dependence 
of foreground contamination on channel frequency in all of 
the WMAP low resolution maps, with the magnitude of 
this excess disappearing in the cleaned, 'foreground reduced' 
maps. 

In the future, a better understanding of the polariza- 
tion foregrounds will allow us to develop more sophisticated 
statistics to precisely evaluate the possibility of foreground 
contamination in cleaned maps. Given the current state of 
understanding, however, we feel that calculation of T> is a 
useful technique for this purpose, as it is one of the sim- 
plest possible statistics which captures the directional signa- 
ture of the polarization foregrounds. It has the advantages 
of being simple and rapidly computable, and can be con- 
sidered complementary to other sta tistics such as the BiPS 
ijBasak. Haiian fc Souradeep 1 2006). Hence it should be one 
of the diagnostic tests ready to assess foreground contam- 



ination in data from the Planck satellite, as well as other 
large-angle polarization experiments. 
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